void diagnose_(double* u, double* ut, double* dx, double* dt);
void diagnose(double* u, double* ut, double dx, double dt){
	diagnose_(u,ut,&dx,&dt);
}
